First principles study of the multiferroics BiFeC>3, Bi 2 FeCrOe, and BiCrO.3: Structure, 

polarization, and magnetic ordering temperature 
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We present results of an ab initio density functional theory study of three bismuth-based multifer- 
roics, BiFeOa, Bi2FeCr06 and BiCrOa. We disuss differences in the crystal and electronic structure 
of the three systems, and we show that the application of the LDA+C/ method is essential to obtain 
realistic structural parameters for Bi2FeCrOe. We calculate the magnetic nearest neighbor coupling 
constants for all three systems and show how Anderson's theory of superexchange can be applied to 
explain the signs and relative magnitudes of these coupling constants. From the coupling constants 
we then obtain a mean-field approximation for the magnetic ordering temperatures. Guided by our 
comparison of these three systems, we discuss the possibilities for designing a multiferroic material 
with large magnetization above room temperature. 



I. INTRODUCTION 



There has been growing recent interest in the magnetic, 
ferroelectric, and piezoelectric properties of bismuth- 
based perovskite-structure oxidesi ^Vi 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 In such 
materials, the formally Bi 3+ ion occupies the perovskite 
A site, and its stereochemically active 6s 2 lone pair 
induces a symmetry-lowering structural distortion that 
can lead to ferroelectricity? 2 * Part of this interest lies 
in increasing concern about the environmental toxicity 
of lead-based piezoelectrics and the resulting search for 
alternative materials In addition, the stereochem- 
ical activity of Bi is being exploited to induce struc- 
tural distortions in magnetic oxides, with the goal of 
forming ferromagnetic ferroelectrics (so-called magne- 
toelectric multiferroics) i^^i^i^ Although a large 
number of potential applications can be envisaged for 
a material which is simultaneously ferromagnetic and 
ferroelectric fi 2 " there are currently no known single-phase 
materials which show large, robust magnetization and 
polarization at room temperature. 

In this work we analyze the magnetic properties of 
three related Bi-based multiferroics in order to under- 
stand the chemical and structural factors which affect the 
magnetic ordering temperature. Our goal is to develop 
guidelines for designing new perovskite-structure ferro- 
electrics with a large macroscopic magnetization above 
room temperature. The focus of our investigation is 
the (111) layered double perovskite Bi2FeCr06, and its 
parent compounds BiFeOa and BiCr03. BiFeOs is the 
most well-studied of the three compounds and has long 
been known to be ferroelectric, with a Curie temperature 
Tc ~ 1100 K, and antiferromagnetic, with a Neel tem- 
perature Tjv ~ 640 Kii&ii It crystallizes in a rhombo- 
hedrally distorted perovskite structure with space group 
R3ci^ Recently, interest in this material has increased 
considerably, due to the report of a large electric polar- 
ization and a non-zero magnetization in epitaxial films 
of BiFeOa, both above room temperature^ The large 
electric polarization, as well as the small magnetiza- 



tion found in 200-400 nm thick BiFeC-3 films, have been 
confirmed by other experiments and by first principles 
calculations^*^ (see also Section lill C|) . whereas the ori- 
gin of the large magnetization observed for film thickness 
< 100 nm is still under debate. 

BiCrOs was recently synthesized and reported to be a 
highly distorted perovskite with C2 symmetry^ a polar 
space group which permits the occurrence of ferroelec- 
tricity. Antiferromagnetic ordering was reported with a 
residual magnetization of 0.02 [Lb per Cr, suggestive of 
weak ferromagnetism. 

Bi2FeCr06 has not yet been realized experimentally, 
but has been predicted theoretically to be ferrimagnetic 
(with a magnetic moment of 2 (is per formula unit) and 
ferroelectric (with a polarization of ~ 80 /jC/cm 2 )^ The 
predicted ground state structure is very similar to the 
R3c structure of BiFeC-3, except that in every second 
(111) layer the Fe 3+ cations are replaced by Cr 3+ , which 
reduces the symmetry to the space group R3. Such a 
-B-site ordered structure could be realized by (111) layer- 
by-layer growth on an appropriate substrate^ 

The remainder of this article is structured as follows: 
First we present the computational details of our cal- 
culations. Then, in Section II I II A-C, we compare the 
calculated ground state structures of the three systems, 
with a particular focus on the influence of the exchange- 
correlation potential used in the density functional for- 
malism on the structural parameters; this is particularly 
important for reproducing the correct physics in these 
strongly-correlated magnetic insulators. We also com- 
pare the calculated electric polarization of the three ma- 
terials. In Section IIII Dl we calculate nearest and next- 
nearest neighbor magnetic couplings and determine the 
magnetic ordering temperatures for the three materials 
using the mean-field approximation. Finally, in Sec- 
tion [W] we discuss the variation of the nearest neighbor 
couplings in these systems and the implications of our 
results for the design of new multiferroics with higher 
magnetic Curie temperatures. 
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II. COMPUTATIONAL DETAILS 

All calculations described in this work were performed 
using the projector augmented wave (PAW) formalism 17 
of density functional theory^iS implemented in the 
Vienna Ab-initio Simulation Package (VASP)^2*2i For 
BiFeC>3 the Bi 5d and the Fe 3p semicore states were 
treated as valence electrons, for Bi2FeCr03 and BiCrC>3 
we did not include semicore states in the valence. The 
local spin density approximation (LDA) as well as the 
LDA+f/ method in the so-called fully localized limit 22 
was used to treat exchange and correlation. Except 
where otherwise noted, the structural parameters were 
fixed to those obtained by full structural optimization us- 
ing U=3 eV and J=0.8 eV (for Bi 2 FeCr0 6 and BiCr0 3 ) 
or J=l eV (for BiFeOs) for the treatment of the transi- 
tion metal d orbitals. Then, the electronic structure and 
magnetic coupling constants were calculated as a func- 
tion of U (and fixed J) in the range U—i eV - 6 eV. For 
simplicity, the same U and J were used on both the Fe 
and Cr sites in Bi2FeCrOg. All our results are well con- 
verged with respect to k-point mesh and the energy cutoff 
for the plane wave expansion. For details regarding the 
structural optimization of BiFeOs see Ref. 0, and Ref. 6 
for Bi2FeCr06- The electric polarization was calculated 
using the Berry phase methodi^Si 

III. RESULTS AND DISCUSSION 

A. Ground state structures and magnetic orderings 

We begin by comparing the calculated structures of 
the three compounds. The ground state structure of 
Bi2FeCr06 was calculated in Ref. |6| by optimizing the 
geometries of a range of starting configurations, obtained 
by freezing in high-symmetry unstable phonon modes of 
BiCrOs 3 and BiA103iiS The energies of structures con- 
strained to the resulting symmetries are listed in Ta- 
bleQJfor two magnetic configurations: ferromagnetic, and 
the ferrimagnetic equivalent of so-called G-type antifer- 
romagnetic ordering, in which all spins within the same 
(111) plane are ferromagnetically aligned, with antipar- 
allel alignment of spins in adjacent (111) planes. Wc 
see that the lowest energy structure of all the combi- 
nations tried has i?3 symmetry (space-group 146), with 
alternating rotations of the oxygen octahedra around the 
[111] direction, combined with relative displacements of 
all ions along [111]. This structure is closely related to 
the R3c structure found experimentally for BiFeOs ^ but 
with an additional symmetry lowering due to the differ- 
ent B cations. From Table [I] one can also see that the 
ferrimagnetic configuration has a lower energy than the 
ferromagnetic case for all structural symmetries. The 
net magnetic moment of the ferrimagnetic case is 2 \xb 
per Fe-Cr pair, which corresponds to a magnetization of 
160.53 emu/cm 3 in the case of the ground state i?3 struc- 
ture. It is also clear from Table [i] that typical energy 



TABLE I: Total energies per formula unit of Bi2FeCr06, for 
different structural symmetries relative to the ground state 
structure for ferromagnetic (FM) and ferrimagnetic (FiM) or- 
dering, respectively. 



space group 


magnetic ordering 


TP 7 /£ . 1 

h [ev/i.u.J 


Pm3m 


FM 


1.722 




FiM 


1.670 


Cmca 


FM 


1.215 




FiM 


1.096 


PA/mnc 


FM 


0.837 




FiM 


0.683 


R3m 


FM 


0.612 




FiM 


0.473 


R3 


FM 


0.176 




FiM 


0.000 



differences between different magnetic configurations are 
significantly smaller than the energy differences due to 
different structural symmetries. Moreover, the struc- 
tural parameters (not shown) obtained for FM and FiM 
ordering within a given structural symmetry differ only 
slightly. 

To exclude the existence of a more complicated mag- 
netic configuration with lower energy, we doubled the 
size of the unit cell along one of the rhombohedral lattice 
vectors and compared the total energies for all possible 
collinear spin configurations within the resulting super- 
cell. The doubled unit cell allows for spin configurations 
with antiparallel alignment of magnetic cations of the 
same type, which would lead to a cancellation of the 
macroscopic magnetization. For these calculations no ad- 
ditional structural relaxations were performed; the struc- 
ture was fixed to that obtained by relaxation within the 
ferrimagnetic i?3 symmetry. Based on the very similar 
structures we obtained above for FM and FiM orderings, 
we do not expect that further structural relaxation would 
alter the relative energies of the different spin configura- 
tions. It was found that the "G- type-like" ferrimagnetic 
configuration described above is indeed the ground state, 
and no cancellation of the macroscopic magnetization oc- 
curs. We did not investigate the effect of spin-orbit cou- 
pling and the possibility of noncollinear spin structures 
or long-wavelength magnetic ordering. 

For BiFeOs both the structure and the magnetic con- 
figuration are well establishedii^ 2 ^ Therefore, we did not 
perform a similar thorough ground state search for this 
system and instead only relaxed the structure within the 
experimentally found R3c symmetry in combination with 
G-type antiferromagnetic ordering. The details of the 
structural relaxation of BiFeOs have already been pre- 
sented in Ref. and all calculated structural parame- 
ters agree well with experimental data. We note that 
the true magnetic structure of BiFeOs shows some slight 
deviations from the ideal G-type ordering, i.e. a long- 
wavelength spiral structure in the bulk2& and a slight 
canting of the magnetic moments in thin films^** both 



TABLE II: Calculated lattice constant a, rhomb ohedral an- 
gle a, volume V, and Wyckoff parameters for BiCrOa, 
Bi2FeCrOe, and BiFeOa. The internal structural parame- 
ters refer to the Wyckoff positions 2a (x,x,x) for the cations 
and 66 (x,y,z) for the anions (In case of the -R3 structure the 
corresponding Wyckoff labels are la and 36). 







BiCrOa 


Bi 2 FeCr0 6 


BiFcOs 


a [A] 




5.43 


5.47 


5.50 


a[°] 




60.37 


60.09 


59.99 


V [A 3 ] 




114.47 


116.86 


117.86 


Bi 


X 


0.000 


0.000/0.503 


0.000 


Cr/Fe 


X 


0.228 


0.226 (Cr)/0.732 (Fe) 


0.228 


O 


X 


0.546 


0.547/0.045 


0.542 




y 


0.952 


0.948/0.450 


0.942 




z 


0.407 


0.405/0.898 


0.368 



caused by spin-orbit coupling. However, these effects are 
rather small and not relevant for the present study. 

In the case of BiCr03 we also constrained the system to 
have rhombohedral i?3c symmetry and G-type antiferro- 
magnetic ordering in our calculations, although recently 
this system was reported to have monoclinic symmetry 5 
This was done in order to be able to systematically com- 
pare the properties of the series of compounds BiFeOa, 
Bi 2 FeCr0 6 , and BiCr0 3 . 

The calculated lattice constants, rhombohedral angles, 
unit cell volumes, and internal structural parameters for 
all three systems are summarized in Table [n] It can be 
seen that the volume and lattice constants increase over 
the series BiCrO^, Bi 2 FeCr06, BiFeOa, due to the larger 
radius of the high-spin Fe 3+ ion compared with Cr 3+ 
(see also Sec. IIII Bp . Likewise the rhombohedral angle 
is reduced. The internal structural parameters, i.e. the 
positions of the atoms within the unit cell are very similar 
in all three systems. This partly reflects the fact that Bi 
is the "active" ion causing the ferroelectric distortion. 



B. Effect of U on the structure of Bi2FeCrOe 

There is an interesting subtlety involved in the struc- 
tural relaxation of the Bi2FeCr06 system, regarding the 
influence of the parameter U on the ground state struc- 
tural parameters. Whereas LDA (U — 0) leads to a 
metallic low-spin state of the Fe 3+ cation (albeit not 
changing the fact that the ground state is ferrimagnetic 
with i?3 symmetry), the use of a moderate U leads to 
an insulating solution containing high-spin Fe 3+ . Due to 
the substantial differences in the ionic radii of low-spin 
versus high-spin Fe 3+ (the Shannon ionic radii are 0.55 
and 0.65 A respectively^), this has a significant effect 
on the lattice parameters. Using LDA for the structural 
relaxation we obtain a lattice constant of a = 5.30 A, 
a rhombohedral angle a — 61.23°, and an equilibrium 
volume that is ~ 7 % smaller than that obtained with 
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FIG. 1: Spin-resolved total and transition metal d densities 
of states (DOS) for U = 0.0 eV (LDA), 1.5 eV and 3.0 eV. 
The total DOS is indicated by the dashed lines, the Fe d DOS 
by shading and the Cr d DOS by thick lines. 



LBA+U and U = 3 eV. 

The evolution of the electronic structure with increas- 
ing U is illustrated in Figure ^ which shows our calcu- 
lated spin-resolved densities of states for U = (the LDA 
result), 1.5 eV, and 3 eV (i?3 symmetry, FiM). To iso- 
late differences in electronic properties from structural 
effects, the same (LDA) structure is used in all cases. 
The orbital-resolved densities of states are obtained by 
projection onto a local angular-momentum basis and in- 
tegrating up to a radius of 1.4 A. In all cases, the spin 
of the Cr ion defines the "up-spin" direction, and its oc- 
cupied {t\ g Y manifold can be seen on the positive y-axis 
just below the Fermi energy. In the LDA solution, the Fe 
t~2 g manifold is completely filled, and the t\ g states are 
partially filled, leading to a finite density of states at the 
Fermi energy. The calculated local magnetic moments 
(also calculated by integrating up to a Wigner-Seitz ra- 
dius of 1.4 A) on the Cr and Fe ions are +2.5 /is and 
—0.6 \ib respectively. These values are slightly reduced 
from the "ideal" values of +3 and —1 expected for Cr 3+ 
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and low-spin Fe 3+ because of hybridization with the oxy- 
gen band and the certain degree of arbitrariness in the 
definition of the Wigner-Seitz spheres. Whereas the po- 
sitions of the Cr states do not change significantly as 
a function of U, the Fe states are strongly affected by 
application of a small U of 1.5 eV; the t\ g states shift 
down in energy and hybridize with the oxygen p band 
(~ 2 eV below the Fermi energy), and the e g states shift 

down into the same energy region as the t 2g states, and 
become partially occupied. On increasing U to 3 eV, 
the Fe e g states become fully occupied, the t\ g states are 
completely empty, and an insulating solution is obtained. 
The calculated local magnetic moments are now +2.5 [Ib 
for the Cr 3+ ion (unchanged from the LDA value) and 
—4.0 hb for the Fe 3+ , again slightly reduced from the 
ideal high-spin value of —5.0 /is- Further increase of U 
does not change the electronic structure qualitatively. 

The above result shows that for Bi 2 FeCr06 the appli- 
cation of the LDA+C/ method induces a major change in 
the electronic structure which, in turn, has a pronounced 
effect on the structural parameters. Traditionally, the 
LDA+C/ method has often been used in combination with 
the atomic sphere approximation (ASA) 28 which does 
not allow structural relaxation. Only recently has the 
LDA+C/ method also been applied to obtain structural 
parame ters, in ge neral with considerable success (see e.g. 
Refs.lfl l29l3(l3ll Our results confirm that, in situations 
where the application of the LDA+C/ method can change 
the spin state of an ion, this latter approach is essential 
for a meaningful result to be obtained. 

Interestingly, if we relax the structural parameters us- 
ing LDA, but fix the total magnetic moment to be 8 
Hb (corresponding to high-spin states for both Cr and 
Fe and a ferromagnetic configuration of all spins), the 
resulting structure is very similar to that obtained by us- 
ing LDA+C/ with U = 3 eV. This shows that the main 
effect of the on-site Coulomb repulsion in L^FeCrOg on 
the structure is the stabilization of the high-spin state 
of the Fe 3+ cation. The larger volume of this high-spin 
state then determines the structure. On the other hand, 
as long as the Fe 3+ cation is in its high-spin state, the 
precise value of U docs not have a pronounced effect on 
the structural parameters. As a consequence, the calcu- 
lated LDA structure for BiFe03 is very similar to the cor- 
responding LDA+U structure, because in both cases the 
transition metal ion adopts the same spin configuration. 7 

C. Calculated polarization 

Table IIIII shows the spontaneous electric polarization 
P s for the three systems, calculated with the Berry-phase 
method (see Section[nJ , as well as an estimate for P s , cal- 
culated by multiplying the ferroelectric displacements, 
Ui, by the formal ionic charges, Zi, and summing up 
over all ions i. It can be seen that, for both methods, 
the polarization increases from BiCr03 to Bi2FeCrOg to 
BiFeOa, reflecting a corresponding increase in ferroelec- 



TABLE III: Spontaneous polarization P a calculated with the 
Berry phase method for all three systems, together with a 
simple estimate, obtained by summing up formal charges Zi 
times displacements if, for all ions i in the unit cell (values 
in /iC/cm 2 ). The ratio between these two quantities is also 
given. 





Ps 


J2i Z i U i 


p s /Y,iZiUi 


BiCr0 3 


67 


46 


1.46 


Bi 2 FeCr0 6 


80 


55 


1.45 


BiFeOs 


95 


64 


1.48 



trie displacements of the ions. The polarization calcu- 
lated from the Berry phase is consistently larger than the 
simple estimate by a factor of ~ 1.5. This represents the 
well-known fact that the effective charges in perovskite 
ferroelectrics are usually significantly larger than the for- 
mal charges. (For example the effective charges in ground 
state i?3c BiFe0 3 are 4.4, 3.5, and -2.5 for Bi, Fe, and O, 
respectively (see Ref. 0) , whereas the corresponding for- 
mal charges are 3.0, 3.0 and —2.0). The nearly constant 
ratio between the simple estimate and the Berry phase 
P s , indicates that the effective charges are comparable in 
all cases. 

The polarizations of all three systems are relatively 
large and comparable to that of PbTiOa (75 /xC/cm 2 )^ 
This suggests that Bi-based perovskite ferroelectrics 
could be promising candidates for alternative "Pb-free" 
ferroelectric and piezoelectric devices^ 



D. Magnetic ordering temperatures and exchange 
couplings 

To estimate the magnetic ordering temperature for 
the newly predicted multiferroic Bi2FeCrOg we deter- 
mine the Heisenberg exchange constants corresponding 
to nearest and next-nearest neighbor magnetic couplings 
by comparing the total energies for different magnetic 
configurations. From these coupling constants we then 
calculate the Curie/Neel temperatures within the mean- 
field approximation (see e.g. Ref. 13311 . For compari- 
son, we perform analogous calculations for BiFeOa and 
BiCr03. Note that the mean- field approximation gives 
an upper limit for the actual ordering temperature, since 
it neglects all fluctuations of the magnetic moments from 
their expectation values. 

To obtain the Heisenberg coupling constants for near- 
est and next nearest neighbor coupling, we double the 
size of the unit cell (so that it contains two formula units 
of Bi 2 FeCrOe) and calculate the total energies for sev- 
eral possible inequivalent collinear magnetic configura- 
tions with different orientations of the magnetic moments 
on different sites. We then map the calculated energies 
onto a simple Heisenberg model, E = —^Y^ijJijSiSj- 
Here Si = ±| on the Cr sites and ±| on the Fe sites, 



t/[eV] 



3 



4 5 
C/[eV] 



FIG. 2: Nearest-neighbor coupling constants, Jij, for BiFe03, 
Bi2FeCr06, and BiCrOs in their ground state structures as 
well as the corresponding coupling constants for Bi2FeCrOs 
and BiCrOs in the ideal cubic perovskite structure as a func- 
tion of U. Dashed lines are for the cubic structures, full lines 
are for the relaxed structures. Circles are for BiCrO^ , squares 
for Bi2FeCr06, and diamonds for BiFe03. 



and the sum is over all nearest and next-nearest neigh- 
bor pairs. All calculated total energies can be reproduced 
accurately using only nearest and next-nearest neighbor 
couplings. 

In Refs. and |3^, the coupling constants for 
Sr 2 FeM0 6 (M=W, Mo, Re) and for BaFei 2 0i 9 were 
found to be rather sensitive to the value of U. This 
reflects the dominant role of the superexchange interac- 
tion in magnetic oxides, in which the antifcrromagnetic 
coupling is proportional to b 2 /U , where b is the transfer 
integral describing hybridization effects^ We therefore 
calculate all coupling constants as a function of U in the 
range 3 eV - 6 eV using the structural parameters shown 
in Table [n] J was set to 1 eV for all systems. To separate 
out effects of changing structure from differences in chem- 
istry, we also calculate the exchange couplings for ideal 
cubic perovskite BiCrO^, BiFeC>3 and Bi2FeCr06, at the 
same average volume of 116.2 A 3 . Figure |21 shows the 
nearest neighbor coupling constants for all three relaxed 
systems and for perfect cubic BiCrOs and Bi2FeCrC>6 as 
a function of U. (For cubic BiFeC>3 some magnetic con- 
figurations are metallic up to U = 6 eV. Since in the 
metallic state different magnetic coupling mechanisms 
are present, we do not include the values for cubic BiFeC>3 
in our systematic comparison.) For both cubic and re- 
laxed structures we obtain a reduction in the strength 
of the antiferromagnetic superexchange as U is increased 
(Figure^; indeed for cubic Bi2FcCrOg the nearest neigh- 
bor coupling even becomes ferromagnetic. 

From the coupling constants Jjy we calculate the mag- 
netic ordering temperature using the mean-field approx- 
imation (see e.g. Ref . l33|l . Figure shows the resulting 
magnetic ordering temperatures for all three systems. 



FIG. 3: Magnetic ordering temperatures for structurally re- 
laxed BiCrOs, Bi2FeCr06, and BiFeC>3 as a function of the 
parameter U. Circles are for BiCrOs, squares for Bi2FeCr06, 
and diamonds for BiFe03. 

Comparing the values of Tjy for BiFeC>3 with the ex- 
perimental value of 640 K shows that for U < 5 eV, 
as expected, the mean-field approximation overestimates 
the magnetic ordering temperature but gives the cor- 
rect order of magnitude. For BiCrO^ and Bi2FeCrOg 
the ordering temperatures are significantly smaller than 
for BiFe03, a fact that can be attributed to the pres- 
ence of strong e g -e g coupling in the latter compound (see 
Section llVjl . Therefore, we expect that the magnetic or- 
dering temperature of the recently predicted multiferroic 
Bi 2 FeCr0 6 will not exceed 100 K. The calculated M(T) 
curve for Bi 2 FeCr06 is shown in Fig. 0] We see that no 
fcrrimagnctic compensation temperature (where the two 
sublattice magnetizations cancel exactly) is observed in 
the temperature range < T < Tq. 

IV. DISCUSSION 

According to Anderson 3 ^ the superexchange interac- 
tion can be separated into two contributions: (i) kinetic 
exchange which is due to the mixing of the ligand-field 
orbitals used to describe the spin quasiparticles; this 
term is proportional to b 2 /U and is always antiferro- 
magnetic, and (ii) potential exchange which represents 
the direct exchange interaction between these ligand- 
field orbitals; this term is always ferromagnetic. Based 
on this separation, Anderson gave the following simple 
guidelines for estimating the signs and relative magni- 
tudes of superexchange interactions (note that the same 
outcomes are reached using the related Goodenough- 
Kanamori rules^£): 

1 . Usually, the kinetic exchange is much stronger than 
the potential exchange and dominates, leading to 
the predominantly antiferromagnetic interactions 
found in magnetic oxides. 
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FIG. 4: Calculated temperature dependence of the magne- 
tization for Bi2FeCrOe (U = 3 eV). The dashed and dotted 
lines show the magnetizations of the Fe and Cr sublattices 
respectively. 



2. The kinetic exchange between e g electrons on dif- 
ferent ions connected by a 180° metal-oxygen-metal 
bond is much stronger than the kinetic exchange 
between correponding t 2g electrons, since the for- 
mer is mediated by dpi bonds whereas the latter is 
mediated by weaker dpi: interactions. 

3. In certain situations the kinetic exchange vanishes 
by symmetry so that the remaining potential ex- 
change leads to a small ferromagnetic coupling. 
This occurs for example for a 90° superexchange 
coupling between ti g electrons and also in the case 
of a 180° coupling between e g and ti g electrons. 

4. Completely filled manifolds, with equal numbers of 
up- and down-spin electrons give no net contribu- 
tion to the superexchange interaction. 

These rules explain, for example, the relative Neel tem- 
peratures of the rock-salt structure transition metal 
monoxides shown in Table ITVl Ni 2+ has electron config- 
uration (tig) 3 , (4 9 ) 3 : ( e J) 2 - Therefore, only the two e g - 
e g interactions contribute to the net superexchange in the 
180° Ni-O-Ni bond (the t\ g - e g and t^ g - e g interactions 
are equal and opposite); these are strongly antiferromag- 
netic, and the Neel temperature is high. For Mn 2+ , with 
electron configuration {t\ g ) 3 , (ej) 2 , there is an additional 
ferromagnetic t<i g - e g direct exchange, which weakens the 
antifcrromagnetic coupling and lowers the Neel temper- 
ature. CoO and FeO are intermediate. 

In BiCrC>3 the balance of the kinetic and potential ex- 
change leads to a moderate antiferromagnetic coupling 
between the ti g electrons on neighboring Cr sites. This 



TABLE IV: Neel temperatures for different transition metal 
oxides. All values except for LaFeOa were taken from Ref. l38fi . 



Material 


T N [K] 


MnO 


118 


Feo.930 


198 


CoO 


289 


NiO 


523 


LaCrOs 


282 


LaFeQ 3 


740 (Ref. 39) 



is stronger in the ideal cubic structure than in the dis- 
torted R3c structure because the perfect 180° Cr-O-Cr 
bond angles in the high-symmetry case result in stronger 
kinetic exchange. 

If half of the Cr 3+ ions in the cubic structure are 
replaced by d 5 Fe 3+ , the antiferromagnetic part of the 
nearest neighbor coupling — now between Fe 3+ (d 5 ) and 
Cr 3+ (d 3 ) — is drastically reduced and in fact the inter- 
action becomes ferromagnetic. This can be explained by 
the additional ferromagnetic interactions between the e g 
electrons of the Fe 3+ and the t2 g electrons of the Cr 3+ 
(see rule 3 above). When U is increased, the remaining 
antifcrromagnetic kinetic exchange between the ti g elec- 
trons is suppressed and the net interaction in the cubic 
structure becomes ferromagnetic. Interestingly, in the 
relaxed structure of R3 Bi 2 FeCr06, the nearest neigh- 
bor coupling between Fe 3+ and Cr 3+ is weakly antiferro- 
magnetic. This difference between the cubic and relaxed 
structures can also be explained within the theory of su- 
perexchange: when the structure is relaxed, the Fe-O-Cr 
bond angle deviates from the perfect 180°, allowing a 
certain degree of mixing between Cr t 2g and Fe e g states. 
This leads to a nonzero kinetic exchange between these 
orbitals and a stronger tendency for antiferromagnetic 
coupling. (Note that similar behavior was discussed in 
Ref. lifl for the La 2 FeCr06 system.) 

Finally, for BiFeC-3 the strong e g -e g coupling due to 
the dpa bonding dominates and results in very strong 
antiferromagnetic nearest neighbor coupling and a high 
magnetic ordering temperature. 

A similar trend in magnetic ordering temperatures 
than discussed here for the BiFe03/Bi 2 FeCr06/BiCr03 
system can also be observed for the closely related 
LaFe0 3 /La 2 FeCr0 6 /LaCr0 6 system (see Table IT?) . 
The ordering temperatures of the end members LaFeOs 
and LaCrOs are comparable to the corresponding Bi 
systems but appear to be slightly larger, which could 
be due to smaller structural distortions in the non- 
ferroelectric La compounds. A Curie temperature of 
375 K and ferromagnetic order were reported for (111)- 
layered La 2 FeCr06, 41 but these results are still under 
debate^S^ From our discussion, ferromagnetic order 
could be expected for this system, although the reported 
Curie temperature seems to be too high to be explained 
by the weak ferromagnetic potential exchange caused by 
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the superexchange interaction. 

From the above discussion it can be seen that high 
magnetic ordering temperatures can be achieved in 
antifcrro- or ferri-magnets by exploiting the strong anti- 
ferromagnetic superexchange between e g electrons. The 
highest ordering temperatures in perovskite structures 
with octahedral crystal field splittings can occur between 
two ions of d 8 electron configuration, although this would 
of course lead to an antiferromagnet with a net cancel- 
lation of magnetic moments. Therefore, we propose that 
the best choice for a ferroelectric with a net spontaneous 
magnetization is to combine d 5 and d 8 ions. Here, pro- 
vided that the crystal field splitting is larger than the 
bandwidth, the material will be insulating, with a net 
magnetic moment of 3 /is per formula unit. While other 
combinations, such as d 8 - d 7 would also give strong su- 
perexchange coupling, the additional partially-filled sub- 
shells reduce the probability of insulating behavior. In 
addition, Jahn- Teller distortions and orbital ordering ef- 
fects can lead to antiferromagnetic coupling in one di- 
rection and ferromagnetic coupling in another direction, 
which makes the simple rules outlined above more diffi- 
cult to apply. 

Finally we discuss some possible combinations of ele- 
ments that will allow d 5 - d 8 coupling in the perovskite 
structure. First we consider retaining Bi 3+ as the A-site 
cation in order to exploit its well-understood stereochem- 
ically active lone pair to induce the ferroelectric struc- 
tural distortion. In this case an average oxidation state of 
3+ is required on the cation B site, which limits the pos- 
sible choices of magnetic ions considerably. Co 4+ /Ni 2+ 
would be a possible candidate for strong d 5 -d 8 coupling, 
but it is questionable if such a high oxidation state can 
be achieved for Co. 

Another possibility is to change the oxidation state 
of the anions, for example by forming oxyfluorides. For 
example Bi2Mn(II)Ni(II)04F2 would have the appropri- 
ate d 5 and d 8 electron configurations, and could give the 
added benefit of an enhanced spontaneous electric polar- 
ization by strategic placement of the F _ anions. How- 
ever, the magnetic coupling via F~ is weaker than via 
2 ~ and oxyfluorides are generally difficult to prepare 
and explosive. 

A better choice would be to change to a different A- 
site cation, which is still lone-pair active, but which has 
a higher oxidation state. For example a 4+ A-site cation 
would retain charge neutrality with divalent d 5 Mn 2+ 
and d 8 Ni 2+ on the B sites. In fact, several perovskite 
systems containing lone-pair active Te 4+ on the A-site 
have been prepared in the past^ 3, with recent reserch fo- 



cussing on Sei_ x Tea;Cu03, as a model system to study 
superexchange interactions^ 4 ^ However these systems 
are usually strongly distorted due to the small size and 
the stereochemically active lone pair of the Te 4+ ion, 
leading to weak kinetic exchange and small magnetic 
ordering temperatures. 43 Po 4+ is also a possibility, but 
since it is radioactive it is less interesting for technological 
applications. 

Moving to the left in the periodic table, Pb 2+ and Tl + 
are also lone pair active cations, and Pb 2+ of course is 
well known as an A-site active cation in ferroelectric per- 
ovskites. Here, however, the average B-site oxidation 
state must be 4+, excluding the possibility of a d 8 config- 
uration. Therefore it is unlikely that room-temperature 
Pb-based perovskite multiferroics will be identified with 
this mechanism. Even disregarding its toxicity, the va- 
lence of Tl + is even more prohibitive. 

In summary, although the spontaneous electric po- 
larization of many perovskite multiferroics persists to 
temperatures far above room temperature, our analy- 
sis suggests that achieving large room temperature mag- 
netization will continue to be challenging. We suggest 
Te2MnNiC>6 as one possible candidate within the per- 
ovskite structure, but it remains to be seen if the struc- 
tural distortions in this system are too large to result in 
strong superexchange coupling. Our analysis indicates 
that other routes to magnetism, such as the weak fer- 
romagnetism believed to occur in BiFe03^ are worthy 
of further pursuit, and that the search for multiferroics 
should be extended to structures other than perovskite, 
in which stronger magnetic interactions might be ob- 
tained. 
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